SUNY-NTG-94/6 



Spectral sum rules and Selberg's integral formula 



Jacobus Verbaarschot 

Department of Physics 
SUNY, Stony Brook, New York 11794 



Abstract 

Using Selberg's integral formula we derive all Leutwyler-Smilga type sum rules for one 
and two flavors, and for each of the three chiral random matrix ensembles. In agreement 
with arguments from effective field theory, all sum rules for Nf = 1 coincide for the three 
ensembles. The connection between spectral correlations and the low-energy effective 
partition function is discussed. 
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1. Introduction 



As has been shown by Leutwyler and Smilga [TJ, chiral perturbation theory of the 
low-energy effective QCD partition function leads to sum rules for the inverse powers of 
the eigenvalues of the Dirac operator in QCD. On the other hand these sum rules can 
also be derived from chiral random matrix theory |2], ||| with only the symmetries of the 
QCD Dirac operator as input. They follow from the microscopic correlations between 
the eigenvalues near zero virtuality, which are conjectured to be universal 0, ||. In || 
we have argued that there are three different universality classes: SU(2)-gauge theory 
with fundamental fermions, SU(N C ), N c > 3 with fundamental fermions and SU(N C ), 
N c > 2 with adjoint fermions. The Dirac operator can be chosen real, is complex or 
can be regrouped into quaternions, respectively. In analogy with the classical random 
matrix ensembles ||, the corresponding chiral random matrix theory will be called the 
chiral orthogonal ensemble (chGOE), the chiral unitary ensemble (chGUE) and the chiral 
symplectic ensemble (chGSE). As was shown in the spectral correlations of chiral 
random matrix theory are different in each of the three cases. Indeed, for two or more 
flavors the effective theory is also different [I], |5|. However, for one flavor, the low-energy 
effective partition function, with proper identification of the parameters, is the same in 
each of the three cases |J , and therefore also the spectral sum rules are the same. This 
leads to the paradox 0] of how an infinite family of sum rules can be the same, whereas 
the spectrum near zero virtuality is different. This issue will be resolved below by an 
explicit calculation of all sum rules. 



In general, spectral correlation functions can be calculated exactly, but the integrals 
for the sum-rules can be evaluated analytically only for the simplest correlation functions 
|3J]. However, there is a more direct way to obtain the same results. The spectral sum 
rules can be expressed in terms of the joint eigenvalue density. The resulting integrals are 
very complicated, but fortunately, exactly these types of integrals have been considered 
before by Selberg H and can easily be generalized with the help of Aomoto's proof H] 



thereof (see |10| for a detailed discussion). Using this technique we are able to derive all 
sum rules for one and two flavors. 
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2. Leutwyler-Smilga spectral sum rules 



For Nf flavors with masses nif (rrif — > 0) the QCD partition function in the sector 
with v zero modes is defined by 

Z? CD (m) = <II II + ™})m}} Sv(Ah (1) 

f=l A n >0 

where the average is over gauge field configurations with v fermionic zero modes weighted 
by the gauge field action S U (A). The product is over all eigenvalues of the Dirac operator. 
For fermions in the adjoint representation, the doubly degenerate eigenvalues count only 
once (Majorana fermions) [|IJ. All sum rules that can be obtained from chiral perturbation 
theory follow by equating the expansion in powers of m 2 of the QCD partition function 
(1) and the low-energy effective partition function. For Nf = 1, the fermion determinant 
can be expanded as 

I[^ 2 n + m 2 )= II A^(l + £™ 2P £ A2 * a V (2) 

A„>0 A n >0 \ p=l m^-^np "l n p J 

Therefore, one finds sum rules for the following quantities 

WlT^-^rij, "i n P I 

where the average, which includes the fermion determinant in the chiral limit, is with 
respect to the QCD partition function. For Nf flavors, it is clear that sum rules for the 
quantities 

\ ^ \ 2N f V^T . 2N f -2 , 2N f -2 7^ / " (^) 

\all A fc different A ai • • • A Q A^ • • • A a • • • A,.„ •••A,. 

* f 



follow from a similar expansion of the product 

n n^+mj). (5) 

A„>0/=1 

Notice that in the thermodynamic limit these sum rules are only sensitive to the spectrum 
near zero virtuality. In other words, they reflect on the microscopic correlations of the 
spectrum, i.e., correlations between eigenvalues near zero on the scale I/V4. 



3. Random matrix theory 



In this paper we calculate the sum rules, not via the chiral perturbation expansion of 
the effective field theory but rather from chiral random matrix theory, which, invoking 
universality arguments 0, [U]], describes the spectrum of the QCD Dirac operator near 
zero virtuality. In this case the partition function is defined by 

Z p , v {m) = J VTPp(T) det ( ^ ^ ) (6) 

where T has the symmetries of the corresponding Dirac operator and the masses are in 
the chiral limit (m./ — > 0). As discussed above, depending on the universality class the 
matrix T is real ((3 = 1, chGOE), complex (j3 = 2, chGUE) or quaternion real (f3 = 4, 
chGSE). In the latter case the square root of the fermion determinant appears in (6). 
The matrix T is a rectangular n x m matrix with \n — m\ = v (for definiteness we take 
m > n), so that the 'Dirac operator' in (6) has exactly v zero modes. The function P(T) 
is chosen gaussian: 

MT) = e W (-^±\£j. (7) 

In this normalization, the chiral condensate equals £ in each of the three random matrix 
ensembles. The total number of modes is N = 2n. The latter quantity is identified with 
the volume of space time. The joint probability density of the eigenvalues of the Dirac 
operator is obtained by transforming to integration variables in which T is diagonal. For 
Nf flavors and topological charge v it is given by 

p,(X u ..., \ n ) = Cf,, n n |A 2 fc - A?f II \l Nf+ ^- 1 exp (-^f £ XI) , (8) 

k,l k \ / k / 

where Cp^ n are normalization constants. 

The expectation values (4) are calculated with respect to the joint probability density 
(8). To evaluate the integrals, it is convenient to introduce new integration variables 

Hk = — g— A fc , (9) 
resulting in the joint spectral density (up to a constant) 

P/3,a(/^l, • • • ,/V)4«l • ■ ■ dfl n = n \fl k - ^lfY[^k eX1 p(-^2^k)dfIl • • • dfl n , (10) 

k<l k k 



where 

2 

Note that for Nf = 1 we have the following remarkable property 

P/3,a(/fyl, • • • , (3Hn) ~ [p/3=l, a (^l, • • • , Pn)]' 3 • (12) 



4. Recursion relations for Selberg's integral 

The evaluation of the spectral sum rules can be reduced to the calculation of moments 

(^•••^W- (13) 

Because all integration variables occur symmetrically, they only depend on the partition- 
ing of the powers and will be denoted by 

H P:a (n p: n p _ 1: --- : n ), (J2 n k = n )- ( 14 ) 

k 

Here, n& is the number of times the power k occurs in (13). By integrating d^pp^fXi, • • • , /i n ) 
we find the recursion relation 

H^ a (n p ,n p - 1 ,---,n ) = (p + a)H^ a (n p - 1, + 1, • • • n ) 

+ ^y/ n^^ny-A (15) 

The last term can be evaluated by exploiting that the the joint probability density is a 
symmetric function with respect to all variables. If we denote / dx 1 dx 2 S{x 1 , x 2 )F(xi, x 2 ) 
by (F)s, we can derive for a symmetric function S(xi, x 2 ) and even values of q the relations 

where we have used that {x\x 2 /(xi — x 2 )) s = 0. For odd values of q a similar relation 
can be obtained 

/ v p-q \ i (9- 1 )/ 2 

(f?^) = ^ (xr fe+i)/2 xr (9+i)/2 ) s + e (ar'ar^- 1 )^ (17) 



In this case the series terminates because 



„k+l k \ I 
X l X 2 \ L l^k^k 



%1 — %2/g 2 



xfx k 2 ) s . (18) 



With the help of these relations all terms in the sum of (15) can be expressed in the 
functions Hp A . As a result we obtain the recursion relation 

(3 p 

Hf3, a {n P i V-i, -• • , n ) = (p + a + -n p ^ 1 + f3^2n p ^ q )H^ a (n p -l,n p ^ 1 + l) 

1 q=2 

Q l(p-l)/2] 

+ — X! n p~(2k+i)H(3, a {n p — 1, n p -k-i + 2, n p _( 2 fc+i) — 1) 
^ fc=i 

p [g/2] 

+ /? J] n p _ g ^ Hp ta (n p - 1, n p _, + 1, n p _ 9+i _i + 1, n p „ 9 - 1), 

(19) 

where the largest integer smaller than x is denoted by [x] , and we used the notation that 
all arguments of H not shown explicitly are the same as in the l.h.s. of the equation. 
In general, this recursion relation is very complicated and an analytical solution seems 
illusive. However for p — 1, 2 a linear recursion relation is obtained which can be worked 
out analytically. For p = 1 the result is |L0[ 



#/3, a ( n i> n o) = (o+ 1 + ^n )#/j, a (ni - l,n ) 



ni — 1 



IJ (a + 1 + I (no + fc))#A«(°> ™o + ni). (20) 

fc=0 ^ 



For p = 2, the result is only slightly more complicated flO 



-ff/3,a(n 2 ,ni,n ) = (a + 2 + + (3n )Hp >a (n 2 - 1, n a + 1, n ) 

ra 2 -l ^ 

= J| (a + 2 + -(rii + /c) + (3n )Hp ;a (0, m + n 2 , n ) 
k=o 1 

n 2 -l o 

= Y[(a + 2 + ^(n 1 + k)+(3n ) 
k=o 1 

ni+n 2 -l a 

x H {a + l + ^{n + l))H^ a {0,0,n + n 1 + n 2 ). 



To obtain the last equality we have used the relation for p — 1. 



(21) 



5. Results for the spectral sum rules 



It is now straightforward to evaluate all sum rules for one and two flavors. First, we 
consider the case Nf = 1. Since all eigenvalues occur symmetrically in (3) we have 

In terms of the integration variables (9), the sum rules become 



n 



S >"UJ W (w-^, ' (2S) 

In order to make contact with Selberg's integral, the new average is with respect to pp, a -i- 
Application of (20) yields 

S P = ( U ) {PnJff ■ (24) 

Remarkably, for Nf = 1 the variable (3 drops out of the equation and the sum rule 
simplifies to 

S P = t — ~~tt [ n ) {nY?y. (25) 

This result constitutes the resolution of the paradox posed in the introduction. In fact, 
we have constructed a one-parameter family of spectra which all give rise to the same 
one-flavor sum rules. Using that the total number of modes is iV = 2n the large n limit 
of this sum rule is given by 

Sp ~ p\(p + v)\2*>' (26) 
which coincides with a result obtained in [jT^j, J] on general grounds involving the anoma- 
lous and chiral structure of QCD. For an arbitrary number of flavors, the large n limit of 
(24) is given by 

{N^Y T(u + l + 2(N f -l)/P) 
p 2>! r(v + p + l + 2(Nf-l)/py 

which again agrees with previous work ||T|. 

The most general sum rule for Nf = 2 is given by 

Spq= \ ^ A 4 • ■ ■ A 4 "'a 2 •■■A 2 )• (28) 

\all \i diflferent A mi /v m p /v ni n n q / 



These sum rules are also well-defined for any number of flavors Nf > 2, and for less than 
two flavors only under certain conditions on v and (3. Therefore, we will evaluate the sum 
(27) for an arbitrary value of Nf, under the assumption that the integral is convergent. 
To evaluate the integrals we change variables as in the case of Nf = 1. Exploiting the 
symmetry of the integration variables we obtain 

-wn—p—q 2 TT n— P 



n fc= i VkRl=n- P - q +iVk) ( BnT?\ q+2p 

OE.,^., (^) ' <29> 

where the total number of terms in (27) is denoted by 

Npq = T7? : TT (30) 

p\q\[n — p — q)\ 

The integrals in (28) follow immediately from the recursion (21). The result can be 

expressed in terms of T-functions: 

o * r r(a + |)r(a+p)r(q + |+p + n) , Y +2 P 

bpq -* pq T(a+l + q + 2p)T(a)T(a + | + n) \" L > ' {61) 

where 

2(N f - 2) 

a = v + l + ^- '-. (32) 

P 

We observe that for p = these sum- rules reduce to the case Nf = 1. For p ^ and 
Nf = 1 the 13— dependence does not drop form (30). This is consistent with the fact that 
in this case the sum rules cannot be obtained from a chiral expansion of the partition 
function. For Nf ^ 1 the result depends on (3. The thermodynamic limit of the sum-rules 
can be obtained with the help of Stirling's formula. This yields 



(iV 2 S 2 )g +2 P r> + f)r(a) 

4<z+2pp!gi r(a + | + g + 2p)r(a+p) 



Two cases that have been considered before in p], namely /3 = 2,p=l,q = and (3 = 2, 
p = 0, q = 2, are reproduced. 



6. Discussion 



In conclusion, we have evaluated all sum rules that follow from a chiral expansion of 
the partition function for one and two flavors. Each sum rule has been calculated for 
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an arbitrary number of flavors and in a sector with a given topological charge. For one 
flavor we have found that all sum rules that can be derived from the effective low-energy 
partition function coincide for each of the three universality classes. This in spite of the 
fact that the spectral correlations are different. This resolves the paradox posed in the 
introduction of this letter. For two or more flavors, both the effective theory and the 
Leutwyler-Smilga sum rules are different for different values of (3. However, we expect 
that also in this case the spectrum is not determined uniquely determined by the effective 
low-energy theory. 

To see the connection between the spectral density near zero virtuality and the low- 
energy effective theory, consider the partition function 

z k v {m, z) = ( n (a* + zr n n ^ + ™?k) • ( 34 ) 

U„>0 /=1A„>0 / Sv ( A j 

The introduction of k replicated flavors allows us to extract information on the spectral 
density. For v = 0, for example, we have the relation 

2np(iz) = lim -^Z^ =0 (m, z) (35) 

To perform the k — > limit of the effective partition function, we first evaluate it for any 
integer value of k and take the limit after analytical continuation. This implies that in 
order to obtain the full spectral density for Nf flavors, we do not only need the effective 
theory for Nf flavors, but for any larger number of flavors (with a different mass) as well. 
Since the effective partition function for an arbitrary number of flavors with equal masses 
is known H, this program may be feasible. 



On the other hand, using techniques described in |Tj], g the spectral density 
near zero virtuality may be related to an effective theory based on graded cosets. At the 
moment it is not clear which approach will be most successful. 
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